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ABSTRACT 

We present an analysis of polarimetric observations at 22 GHz of the water vapor masers in 
NGC 4258 obtained with the VLA and the GBT. We do not detect any circular polarization in 
the spectrum indicative of Zeeman-induced splitting of the maser lines of water, a non-paramagnetic 
molecule. We have improved the l-a upper limit estimate on the toroidal component of the magnetic 
field in the circumnuclear disk of NGC 4258 at a radius of 0.2 pc from 300 mG to 90 mG. We have 
developed a new method for the analysis of spectra with blended features and derive a l-a upper 
limit of 30 mG on the radial component of the magnetic field at a radius of 0.14 pc. Assuming 
thermal and magnetic pressure balance, we estimate an upper limit on the mass accretion rate of 
~ 10~ 37 M Q yr~ 1 for a total magnetic field of less than 130 mG. We discuss the ramifications of our 
results on current maser models proposed to explain the observed maser emission structure and the 
consequences for current accretion theories. We find from our magnetic field limits that the thin-disk 
model and the jet-disk model are better candidates for accounting for the extremely low-luminosity 
nature of NGC 4258, than models that include advection-dominated accretion flows. 

Subject headings: Accretion, Accretion Disks — Galaxies: individual (NGC 4258) — 

Galaxies: Magnetic Fields - Galaxies: Nuclei Masers - Polarization 



1. INTRODUCTION 

Magnetic fields appear to be important for many as- 
pects of material-inflow and -outflow in active galac- 
tic nuclei (AGN), which are powered by accretion 
onto the central super-massive black holes (SMBHs). 
Magnetic fields are thought to provide outward an- 
gular momentum transport in the accretion process 
via the magneto-rotational-instability (MRI), with ac- 
compa nied dissipative heating, allowing furt her mass- 
inflow (H2Eni^!l^Ii3 ElM ISalEuatinQS)- Strong 
poloidal magnetic fields may thread centrifugal winds, 
which could carry away angular momentum, while 
confining and collimating la rge-scale outflows and jets 
seen in radio-loud AGN s l|Blandford & Pavnei Il982t 
iKrasnopolskv et alJll999j) . Furthermore, broad-line re- 
gion (BL R) clouds m ay be constrained by magnetic 



winds (e.g.. 
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NGC 4258 is a very good laboratory in which to 
attempt to measure magnetic fields, since the mag- 
netic field in an accretion disk could be directly mea- 
sured very close to the central black hole. NGC 4258 
(M 106) is a low-luminosity, weakly active Seyfert II 
galax y (or of Seyfert type 1.9, as defined bv IHo et al.1 
1997). It possesses an obscured X-ray center with 
an absorbing column of Nh ~ 10 23 cm~ 2 and an 
absorption corrected 2 — 10 ke V luminosity of ~ 
10 41 erg s -1 ijFruscione et al.ll2005[) . Nuclear continuum 
and narrow emiss ion lines have been detected in polar - 
ized optical light {Wilkes et al.lll995tlB~arth et al.l Fl999L 

The discovery of a nu clear megamaser in a survey 
by iClaussen et all l)1984|) in NGC 4258 and the subse- 
que nt radio single-dish and interferometry ima ging stud- 
ies l|Greenhill et alJl!995t iMivoshi et all 1 19951) have im- 
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proved our knowledge of the nuclear region and its kine- 
matics on sub-parsec radii. The spatial distribution and 
line-of-sight velocities of the masers delineate the molec- 
ular portion of a thin, Keplerian (better to 1%) accre- 
tion disk that is ne arly edge-on (tipped down by ~ 8°) 
and s lightly warped l)Moran et alJll99"5llHerrnstein et alJ 
1996). The maser lines fall into three groups: a dense 
cluster near the systemic velocity of the galaxy, which 
lies close to the line-of-sight to the black hole (here- 
after referred to as "systemic maser lines"), bracketed 
by two groups of less densely packed lines offset by ± 
(800—1000) km s _1 . These high- velocity masers have 
very s mall accelerations compared with the systemic 
masers l|Bragg et al.ll2000|) and are distributed over dif- 
ferent radii close to a line perpendicular to the line of 
sight, the so-called "midline" (Fig. [TJ; hence their line- 
of-sight velocities are close to the Keplerian orbital ve- 
locities. The inner radius of the maser annulus is ap- 
proximately 0.14 pc and its outer radius approximately 
0.28 pc (for a distance of D = 7.2 Mpc, Herrnstein et 
al. 1999; 1" corresponds to 35 pc at that distance) from 
the central black hole with A/ B h = (3.9 ± 0.4) x 10 7 M Q 
ijMivoshi et alJl!995|) . There is a prominent north-south 
twisted jet on kiloparsec scales observed in the radio, in 
Ha, and in X-rays (e.g., I?^l2000t IWilsou et a,lJl20f)T1 
and references therein). It appears to be an extension 
of the radio nuclear je t seen at sub-parsec di stances from 
the dynamical center ijHerrnstein et al.ll997|) and to have 
an axis that is well aligned with the rotation axis of the 
maser accretion disk. We note that the angular momen- 
tum vectors are quite misaligned at kiloparsec and sub- 
parsec distances from the central engine: the edge-on 
accretion disk, as traced by the masers, makes an an- 
gle of ~ 120° to the galactic disk of the host galaxy, a 
SAB(s)bc spiral (iTullvl 119881) . i.e., the maser accretion 
disk is counter-rotating with respect to the larger galac- 
tic structure. 
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Fig. 1. — Correspondence between maser features in velocity at 22.2 GHz and their loca tion in the disk. M iddle: The warped accretion 
disk, defined by the kinematics of the maser spots, is shown as a wire grid (adapted from Moran ct al. 1999), while the black dot marks 
the dynamical cente r of the system (the pos ition of the SMBH). The gray-scale elongation in the center shows the 22 GHz synchrotron 
continuum emission IHerrnstein et al.lll998alh Top, bottom: The GBT spectra were taken on 2003 October 23 (see §E|for details). 



There have been estimates made of magnetic 
field strengths at kiloparsec radii in NGC 4258. 
IKrause fc Lohrl lEffll obtained VLA 8.44 GHz (A3. 6 cm) 
continuum observations of the jets of NGC 4258 and de- 
tected linear polarization at a level of 35% — 65% along 
the northern jets whose interpretation in the context of 
synchrotron emission mechanisms suggests the presence 
of an ordered magnetic field. Assuming energy equiparti- 
tion between the magnetic field and the relativistic parti- 
cles in the northern jet, and an electron-proton (electron- 
positron) composition, they obtained an average B field 
of ~310 [iG (^90 [iG). These values reflect the poloidal 
component of the global magnetic field at a distance of 
^1 — 4 kpc along the jet. 

By searching for Zeeman-induced circular polarization 
in the maser lines, the strength and direction of the mag- 
netic field parallel to the line of sight across the maser 
emission region can be directly determined. The effort to 
scrutinize a well-isolated, 1.2 km s _1 w ide maser line in 
the red shifted group at 1306 km s _1 bv IHerrnstein et all 
( 1998b}) yielded no detection of polarization in the VLA 
data taken in 1995. The resulting limit to the circular 
polarization of 0.5% corresponds to an upper bound of 
approximately 300 mG for the toroidal component of the 
magnetic field. The Zeeman technique constitutes the 
most direct way of obtaining the magnitude and direction 
of the magnetic field. Magnetic field estimates have bee n 
derived from other techniques (e.g., Uones e£a0|2003), 



but they are relatively indirect and model dependent, 
relying on assumptions of equipartition between ther- 
mal and magnetic energy, and adopted values of electron 
density and temperature at various radii. The initial at- 
tempt to measure p olarization in NGC 4258 was made by 
iDeeuchi et alJ l)1995|) who reported an upper limit on lin- 
ear polarization of 20% for the high- velocity water maser 
components and of 15% on the systemic components. 

The format of this paper is as follows: In §[3J we briefly 
outline the theoretical background for understanding the 
Zeeman effect in the water molecule, and describe the 
technical details of how to measure the Zeeman effect. 
We present and discuss two data sets in § |3| and report 
the measured upper limits on the magnetic field for the 
various maser lines in § We discuss the interpretation 
of the Zeeman observations in § El an d give conclusions 
with §|S| 

2. THE ZEEMAN EFFECT IN WATER 

The maser emission at 22.235 GHz involves the ro- 
tational energy states Jk+k- — ^23 and Jk+k- = 6i6 
of ortho-water, which lie ^600 K above the ground 
state. This rotational transition consists of six hy- 
perfine components, but it is generally believed that 
only the strongest three hyperfine components (F — 
7 — 6, 6 — 5,5 — 4) contribute significantly to the maser 
emission llDeguchi fc Watsonlll986t iNedoluha fc Watsonl 
Il99lt iVlemmings et al.ll2001l hereafter I VP if with veloc- 
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ity separations of 0.45 and 0.58 km s 1 between them 
(those separations are marginally within a typical mascr 
line width of ~ 0.8-1 km s" 1 ). 

In the presence of an external magnetic field, each hy- 
perhne component splits into three groups of lines with 
characteristic polarization properties: the tt transitions 
(AMj? = 0) produce radiation linearly polarized along 
the external B field, with negligible frequency displace- 
ments from the parent-hyperfine transition. The 
transitions (AM p = ±1) generate emission whose E field 
vector rotates in a circle in the plane perpendicular to the 
external B field direction and are symmetrically offset 
from the parent-hyperfine in frequency by ±A.i>z (we fol- 
low the notation used in lVOlt see their Fig. 1 showing the 
theoretical Zeeman pattern). Since the water molecule is 
non-paramagnetic, the splitting is very small (Avz ~ 1 
Hz mG" 1 or 10 -5 km s _1 mG" 1 compared to, e.g. up to 
3 kHz mG -1 or 0.6 km s" 1 mG" 1 for the OH molecule 
at 18 cm) and much less than the thermal line width 
of the maser features for expected sub-Gauss B fields 
(Av z ~ 10" 4 - 10" 3 A^ L ). For a B field along the line 
of sight, the components appear in left-circular po- 
larization (LCP) and right-circular polarization (RCP), 
respectively, and their slight frequency separation can 
be detected as an "S" shaped profile in the Stokes V 
[= (LCP — RCP)/2 ] spectrum for a simple Gaussian 
line profile (Fig. [2}. This S profile is proportional to 
the derivative of the Stokes I [= (LCP + RCP)/2 1 spec- 
trum in the sma ll splitting case ilFiebig k Oustenlll989l 
hereafter IFG89I: iHerrnstein et alJll998bt IVOH see also 
Appendix |A"|) . 

By modeling the theoretical H2O Zeeman pattern and 
convolving it with Gaussian profiles of differing line 
widths, IFG89I constructed a set of synthetic Stokes V 
spectra and found the following relationship between the 
applied B field and the resulting circular polarization: 

Vmax _ ^F-F'-B|| ^ 
/max Av L 

where V max is the maximum of the S curve, / ma x is the 
maximum in Stokes /, Bu is the magnetic field parallel 
to the line-of-sight, Bn — Bcos9 (where 9 is the angle be- 
tween the B direction and the line of sight), and Af-f' 
is a coefficient that depends on which hyperfine compo- 
nent is assumed. In the case where the hyperfine lines 
appear individually, Aj^q = 0.0134 km s" 1 G" 1 for the 
F = 7 — 6 hyperfine that possesses the largest transition 
probability; A 6 _ 5 = 0.008 km s" 1 G" 1 for F = 6- 5 and 
A 5 _ 4 = 0.001 km s- 1 G- 1 for F = 5 - 4. 

Equation Q does not account for the non-thermal 
nature of maser emission, which includes inverted, 
non-thermal, populati on distributions and amplifica- 
tionof e mission. iNedol uha fc Watsonl Il992l (hereafter 
NW92T). IWatson fc Wiebd 120011 a"nd T more recently, 
Vlemmings et al. [|2C)02 (hereafter lV02fl . conducted a se- 
ries of studies in which they solved the radiative trans- 
fer and rate equations numerically for water masers in 
a magnetic field, extending the pione ering analytic work 
of IGoldreich. Keelev. fc Kwanl IU971) . These numerical 
simulations account for saturation effects and the pres- 
ence of the three strongest hyperfine components and 
their magnetic substates. Their major findings were 
that the "thermal" equation is an excellent approx- 



imation for the case of maser polarization induced by 
a magnetic field with an effective coefficient Ap-F' of 
0.020 km s _1 G _1 (appropriate for the presence of the 
three strongest hyperfines). The estimate of the B field 
via equation Q is accurate to within a factor of 2, as 
long as the water maser is not hig hly saturated, i.e., has 
a linewidth < 0.8 km s~ 1 l|NW92fl . In those numerical 
simulations, maser propagation was modeled i n one di- 
mension. However, initial calculations in 2D bv lV02l sug- 
gest that the resulting Zeeman S curves might be slightly 
narrower and ha ve a distorte d profile, than dedyced from 
ID calculations. IGravl (|2003) did a study of maser polar- 
ization models in which he compared his general mod- 
els to the ones put forth by Watson an d collaborators 
l|Watsonll994HNedoluha k Watsonll9"94]) and found that 
his equations can be reduced to the ones by Watson, for 
the limiting case of small Zeeman splitting. 

In summary, we adopt Af_f' = 0.020 km s _1 G _1 for 
the S'-curve model, where the Stokes V profile is pro- 
portional to the derivative of Stokes /, appropriate for 
the presence of the three strongest water hyperfines in a 
masing environment. This Zeeman splitting coefficient is 
larger by a factor of 1.5 than for the LTE model where 
solely the hyperfine component with the largest transi- 
tion probability under LTE is present. Follo wing other 
Zeeman experime nts in the radio regime fe.g.. lGoodmanl 
1 1 9 8 9llHeilesl2 00 1|) . we adopt the IEEE convention for cir- 
cular polarization, where a magnetic field pointing away 
from the observer is positive and produces a LCP com- 
ponent at a lower velocity than the RCP component. 

Several water maser polarization experiments have 
been conducted and interpreted in the framework of 
equation Q, with line-of-sight magnetic fields inferred 
on the order of 10 — 50 mG in Galac tic star-forming 
regions l)FG89t iSarma et al.ll2001l |2002|) . These results 
agree with magnetic field strengths expected for a wide 
range of probed densities in molecular clouds, which are 
subject to flux freezing and ambipolar diffusion, thus pro- 
viding corroborating evidence for the cogency of the Zee- 
man interpretation. Furthermore, line-of-sight magnetic 
fields have been inferred on the order of 50 — 300 mG in 
circumstellar envelopes of late- type stars l)V02|) . These 
values are consistent with B field estimates based on OH 
maser observations for a reasonable radial scaling of the 
B field. 

There are two types of maser lines in the spectrum of 
NGC 4258, which require slightly different treatments to 
measure the Zeeman effect. For the well-isolated red- 
shifted and blueshifted high- velocity maser lines, the 
standard S-curve analysis method as outlined above 
works well. For the systemic features, which are heavily 
blended in velocity space, we developed a variant of the 
S'-curve method. Our new technique is more appropriate 
since it does not require the a priori knowledge of the ve- 
locity components. In this procedure, we cross correlate 
the LCP and RCP spectra to determine the net velocity 
offset (see Appendix 151 for details ) . This cross correla- 
tion method (see also IHeiled IT 9 8 81 for a slightly different 
treatment in the case of observations of blended thermal 
H I features) is based on the assumption that the B field 
strength and direction remain constant over the emis- 
sion region, i.e., that the Zeeman- induced velocity shift 
is constant across the spectrum. 

From simple statistical arguments, supported by 
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Fig. 2. — Theoretical Stokes V spectra for the three strongest hyperfine components F = 7 — 6, 6 — 5 and 5 — 4, under the assumption 
of LTE at the predicted hyperfine rest velocities (centered at the hyperfine rest velocity of F = 6 — 5). The B field is 50 mG and the line 
width 0.8 km s- 1 (FWHM). 



Monte Carlo simulations C Appendix LB.2f> . the rms sensi- 
tivity to the LOS B field (as,,) is given by the following 
relationship (which applies for both the S'-curve and the 
cross correlation method) , 

_ 1 Av L 1 

B " 2A F _ F< SNR y^/V ' 1 ' 

where Av L is the FWHM line width, A F _ F , the Zee- 
man splitting coefficient, SNR the signal-to-noise ratio 
for the specific line (in Stokes I) and N the number of 
lines over which the fit is performed. This result is ap- 
propriate if the N lines are not heavily blended. For 
differing line strengths and line widths of the spectrally 
blended lines, as found in the cross correlation technique, 
the Avl/ (SNR V~N) factor must be replaced by an ap- 
propriate average. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. VLA Data 

Data were taken on 2001, December 16 in spectral line 
and full polarization mode with the Very Large Array 
(VLA) of the National Radio Astronomy Observatory 
(NRAO) 1 in its D configuration for approximately 12 
hours at K Band (22.2 GHz). The resulting synthe- 
sized half-power beamwidth is 2.8 arcsec. In this beam, 

1 The National Radio Astronomy Observatory is operated by As- 
sociated Universities, Inc., under the cooperative agreement with 
the National Science Foundation. 



the maser spots (with intrinsic sizes of less than 100 
/larcsec and inter-maser-separation s of several milliarc- 
sec as shown by V LBI observations: iMivoshi et al]ll995t 
iMoran et all [1995) appear as one point source and are 
not resolved; therefore, beam effects such as beam squint 
and beam squash do not affect our data. The segrega- 
tion among the various maser spots occurs in velocity 
space, since the maser spots possess different orbital ve- 
locities. The observations of an isolated high-velocity 
feature at 1269 km s _1 had to be interlaced with those 
of the systemic features due to the limited bandwidth 
of the VLA. The calibrators were observed separately 
for the two different frequency bands. The observations 
were performed with a bandwidth of 1.56 MHz with 128 
channels, yielding a channel spacing of 12.2 kHz (or 0.16 
km s _1 ). With uniform weighting the spectral resolution 
is 0.19 km s" 1 , which corresponds to ~ 1/5 of the high- 
velocity maser line width. The first three hours of the 
observations experienced snow fall, with skies gradually 
clearing thereafter. 

The data were edited, reduced, and calibrated with 
standard spectral line techniques in AIPS 2 , including 
bandpass, amplitude and phase calibration. Amplitude 
calibration was performed with respect to 3C 147, a VLA 
calibrator with known flux density at 22 GHz, phase cali- 
brations relative to 1146+399, and bandpass calibrations 
using observations of 0854-201. The system tempera- 

2 Astronomical Image Processing System of NRAO. 
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Fig. 3. — GBT total power data with a bandwidth of 200 MHz. Note the flat baseline after removing a running mean that excluded the 
line-displaying regions. Inserts show zooms on the red- and blue-shifted high-velocity lines, respectively. The rms noise in each channel 
is 2.9 mjy. The feature at 1647 km s _1 is a highly significant new discovery. It is the most redshifted high-velocity feature discovered in 
NGC 4258, and it extends the inner boundary of the maser disk to 0.12 pc, based on the Keplerian rotation curve. 



tures of the K-band receivers in nine of the antennas were 
in the range of 100 — 130 K, whereas the new K-band re- 
ceivers, residing in 17 out of the 26 operational antennas 
at the date of observations, had system temperatures of 
60 — 80 K. The amplitude calibration is estimated to be 
accurate to within 15%. Strong internal radio-frequency- 
interference was present in the calibrator data taken with 
the old receivers at the frequency corresponding to that 
of the systemic maser emission and had to be removed 
from the data set. 

We performed linear polarization calibration in two 
steps: first we used the observations of the strong source 
0854+201 over a range of parallactic angles to correct 
for the imperfect feed responses (as expressed by the so- 
called "leakage terms" , which are on the order of a few 
% for the VLA K band feeds). Secondly, we determined 
the phase offset between LCP and RCP by observing 
3C286, a polarized source with a known and invariant 
polarization angle. For circular polarization calibration, 
an additional step involved correcting for the differences 
in gain for LCP and RCP and was performed during the 
analysis of the polarization data (see § . 

3.2. GBT Data 

We took spectral line data in RCP and LCP on 2003 
October 23 with the Robert C. Byrd Green Bank Tele- 
scope (GBT) of the NRAO. With an off-axis feed arm, 
the GBT delivered an unblocked aperture and low side- 
lobes with a FWHM beamwidth of 33" at 22.2 GHz. 



The spectrometer had a bandwidth of 200 MHz config- 
ured to provide a channel spacing of 0.164 km s^ 1 with 
its 16384 channels over 200 MHz (2700 km s" 1 ). The 
K-band receiver for the upper frequency range (22—26.5 
GHz) had two beams, each beam producing two polar- 
ization channels. The data were taken in beam-switched 
mode with beam throws of the order of 4', the angu- 
lar distance between the pair of beams. The cycle time 
of two minutes was sufficient to compensate for atmo- 
spheric fluctuations. In this fashion, one beam of the 
telescope was always pointing at the source. A total of 
eight hours of data were taken, including pointing and 
focusing observations after every hour of source observa- 
tion. Applying the local pointing corrections improved 
the pointing uncertainty to < 2". The effective on-source 
time was 348 minutes and the zenith opacity ~0.07. 

The gain curve, which describes the efficiency of the 
antenna as a function of elevation, had been determined 
from observations of the radio source NGC 7027 over 
a range of elevations. The gain curve with the active 
surface-control enabled indicates a peak antenna effi- 
ciency of 1) = 0.6 at an elevation of ~ 30° and stays 
comparatively flat (77 ~ 0.5 — 0.6) over a wide range of 
elevations (20°— 80°). The rms deviation of the gain mea- 
surements around the best-fit polynomial model, which 
is applied to the data for gain correction, was about 10%. 
Amplitude calibration was performed via injection of a 
signal from a noise diode every 0.25 sec. Data reduction 



6 Modjaz et al. 



TABLE 1 

Journal of H2O Maser Observations 



Parameter 


VLA Datat 


GBT Datat t 




Set I 


Set II 




Date 


2001 Dec 16 


2001 Dec 16 


2003 Oct 23 


Radio LSR center velocity 3, [km s _1 ] 


495 


1268 


550 


Total Bandwidth [km s" 1 ] 


21 


21 


2700 


Channel spacing [km s -1 ] 


0.16 


0.16 


0.16 


Integration time b [min] 


241 


231 


348 


1 <r noise in Stokes I c [mjy] 


2.4 d 


2.4 


2.9 



tRR, LL, RL, LR observed. 
t+RR and LL observed. 
a V hc i = VisR - 29.169 kms- 1 . 

b Effective time spent on source, not including time spent on calibrators. 
c Noise per channel, corrected for variable antenna gain and referenced to outside 
the atmosphere. 
d Derived from line-free intervals at high velocities. 



was performed with a customized IDL procedure. 

The residual baseline variations in the spectra were 
eliminated by subtracting a running mean computed over 
intervals of 256 channels (41 km s" 1 ) that excluded the 
line-bearing regions; this procedure yielded flat baselines 
in the final spectrum (Fig. |3J) . The system temperatures 
for the GBT data ranged from 35 — 50 K and quoted 
flux densities are accurate to 10%. The largest source of 
uncertainty is the inaccurate knowledge of the strength 
of the noise diodes over the large bandwidth of the GBT 
data. We note that the uncertainty in the absolute in- 
tensity scale has little impact onto the reported results of 
the magnetic field upper limits: only a differential mea- 
surement, the SNR, determines the upper limit (eq. [5]). 

Table ^ lists a summary of the specifications of both 
data sets. The velocities in this paper follow the radio 
definition of Doppler shift {/S.v/v — —v/c) and are com- 
puted with respect to the local standard of rest (LSR). 

4. RESULTS AND INTERPRETATION 

For circular polarization data, instrumental gain dif- 
ferences between LCP and RCP cause a portion of the 
Stokes / profile to appear in the Stokes V spectrum. This 
effect is mitigated by simultaneously fitting a replica of 
Stokes I to the Stokes V spectrum in addition to a scaled 
derivative of Stokes J in the least-squa res sense (e.g., 
ITroland fc HefiellLosl ISault et alJlToM) : 

V(v)=Ci^P-+c 2 I(v). (3) 

Here, C2 denotes the instrumental polarization term. The 
factor ci is a measure of the Zeeman splitting and corre- 
sponds to Af-f'-E>cos#/\/2 f Appendix lA"|l . This calibra- 
tion procedure is equivalent to the assumption that equal 
flux should be present in both the LCP and RCP spectra, 
i.e., that the circular polarization spectra are intrinsically 
anti-symmetric. However, theoretical maser polarization 
simulations (see § produce V spectra that are asym- 
metric due to the presence of hyperfine structure. Efforts 
to verify the intrinsic shape of the V spectrum observa- 
tionally are hampered by the the fact that it has been 
impossible to determine the instrumental circular polar- 
ization independently well enough to the requisite degree 



of P v = V/I ~ 0.1% (IZhao. Goss. fc Diamond! fl993h 
In order for the derived magnetic fields strength to 
be counted as a detection, we require the field strength 
to be greater than 3 times the associated uncertainty, 
<JB\ \ ■ With this criterion, we did not make any Zeeman- 
induced circular polarization detection. 

4.1. High-velocity Features 

FigureHshows the Stokes /[= (LCP+RCP)/2], i.e., to- 
tal power, VLA spectrum of the redshifted high-velocity 
feature at 1269 km s _1 with a thermal rms noise of about 
2.4 mjy (in a 0.16 km s _1 resolution interval; see Ta- 
ble 0). The high- velocity line was best fit by a two- 
component Gaussian, such that the main line is centered 
at 1268.6 km s^and has a FWHM line width of Av L 
— 1.01 km s _1 . A formal minimum y 2 fit of the Stokes 
V[= (LCP — RCP)/2] spectrum to equation J2J indicated 
a 1-er upper limit to any intrinsic circular polarization of 
Pv ^ 0.08% in the 1269 km s _1 feature. For a maser 
Zeeman splitting coefficient of v4p_F' = 0.020 Hz mG -1 , 
cjbh w 40 mG (1 a). This value is the same, as expected, 
as the value derived fro m eq uation J5J and Monte Carlo 
Simulations (Appendix IB.2J) . From a conservative es- 
timate for the upper limit on the magnetic field along 
the linc-of-sight, I B\\ m |=| Bjj utput | +a B]] + a a we 

obtained B\\ < 90 mG (Table EJ), where cr<r B|| stands 
for the uncertainty in o~B n and is estimated from Monte 
Carlo simulations (§ IB. 2(1 . This feature is located at 
a radial distance of 0.27 pc from the SMBH (70,000 
in units of Schwarzschild radii, Rs = 2GA/bh/c 2 , for 
Mbh = 3.9 x 10 7 M Q ). Following convention, the middle 
panel of Figure 01 shows the Stokes V spectrum intrinsic 
to the source, i.e., after removing a residual replica of 
Stokes /. Superimposed are the 1-a Zeeman S curves for 
Bjj m = ±90 mG. The bottom panel of Fi gure 0]shows the 

total linear polarization [Stokes (Q 2 + [7 2 ) 1 / 2 ] spectrum: 
no fractional linear polarization [Stokes (Q 2 + U 2 ) 1 / 2 / 1] 
at a level of 0.1% (1 a) is detected. 

For the GBT dataset, the combination of low sys- 
tem temperatures and a new spectrometer with a large 
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bandwidth promised higher sensitivity for Zeeman de- 
tection. However, the variable nature of masers caused 
a decrease in the SNR during our observations and no 
Zeeman-induced circular polarization could be detected 
in any of the numerous features. The S'-curve analy- 
sis yielded an upper limit for the strongest high-velocity 
feature at 1303.3 km s" 1 of B|| <90 mG, similar to 
the VLA experiment (Table Figure [5] shows the 
Stokes / spectrum (top) with a thermal rms noise of 2.9 
mJy (in a 0.16 km s _1 resolution interval). This fea- 
ture is best-fit by a one-component Gaussian model of 
width 0.79 km s" 1 (FWHM), centered at 1303.3 km s" 1 . 
The instrumentation-corrected Stokes V spectrum is dis- 
played in the bottom panel of Figure[5]with the overlayed 
1-cr S curves for Bj| m = ±90 mG. This feature is located 
at a distance of 0.25 pc from the central black hole, which 
corresponds to 62,000 Rs- Linear polarization data could 
not be obtained with the GBT spectrometer at the time 
of observations. 

In this case, B|| corresponds to the toroidal magnetic 
field and our limits are significant on the assumption 
that B|| is not tangled over the maser emission region 
of < 0.003 pc (or < 10 16 cm) <|Moran et al.lfl995^ . If the 
22.2 GHz water vapor line consists solely of the F = 7 — 6 
hyperfine component, the B field estimates have to be 
multiplied by a factor of 1.5. This result is an improve- 
ment of the upper limit of B| < 300 mG derived by 
iHerrnstein et al.l l)1998b|) on a feature at 1306 km s _1 
using a similar method. 



Hydro dynamic and MHD turbulence are expected in 
accretion disks. Assuming a Kolmogorov power spec- 
trum for turbulent magnetic fields in a circumnuclear 
disk, iWatson k, Wiebel l)2001|) produced realistic simula- 
tions of total-power water maser spectra and found that 
the consequences for circular polarization were not seri- 
ous. In the case of turbulent fields, they suggested that 
equation Q underestimates the magnitude of the mag- 
netic field generally by a factor of about 3. In their calcu- 
lations, they assumed that the turbulence length for the 
radial scale is the same as the scale height of the disk, H, 
as they assumed isotropic turbulence where the longest 
wavelength is set by the height. However, C. Gammic 
(2004, private communication) estimated that most of 
the power may be contained on scales of the order 40-ff. 
Since the length of the masers and H/R are estimated 
to be <0.003 pc and <0.0025, respectively l)Moran et alJ 
1995), the turbulence length according to Gammie is 
larger than the maser length, by a factor of 4—10. For 
turbulence lengths that are considerably larger than the 
maser amplification length, few magnetic field reversals 
are expected across the maser length. Thus in summary, 
we conclude that it is unlikely that the non-detection 
of Zeeman-induced circular polarization is due to a sce- 
nario where turbulent magnetic fields quench the maser 
polarization signal due to magnetic field reversals. 

4.2. Systemic Features 

For the cross correlation method, we performed an 
equivalent gain calibration procedure by requiring the 
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integrated line fluxes of LCP and RCP to be equal. Sim- 



ulations and application of both calibration methods to 
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TABLE 2 

H2O Maser Parameters and Zeeman Interpretation 



Data Set Line center a FWHM I D Fitted Range B|| b cr S|| b a aB c Bjj m (lcr) d R 
[kms- 1 ] [kms- 1 ] [mjy] [km s" 1 ] [mG] [mG] [mG] [mG] [pc] 



VLA e 1268.6 1.01 1795 1266.5-1270.0 -47 38 5 90 0.27 s 

GBT e 1303.3 0.79 1212 1302.0-1304.5 40 44 6 90 0.25 g 

VLA f 499.1 ... 4553 486.0-504.0 29 16 4 49 0.14 h 

GBT f 492.3 ... 4102 472.0-507.0 7 19 4 30 0.14 h 



a Radio LSR velocity. 

Values for _B|| and gg^ have been obtained from the Zeeman equation with Ap_ F / = 0.020 

kms" 1 G" 1 INedoluha fc Watsonl [19921) . which presumes that the 22.2 GHz water maser line is a merger 
of the three strongest hyperfme lines. The values for crg^ agree with values derived from Monte Carlo simula- 
tions with characteristic noise and data parameters (see 5 IB.2I . 
c Values for a aB were obtained from Monte Carlo simulations with synthetic data that display the same data 

and noise characteristics as the corresponding datasets (see S IB. 21 . 
d Obtained from a conservative estimate for Bjj m = | £j° ut P ut | _|_ ag^ ^^^b^ ■ ^ the ^2.2 GHz water vapor 

line consists solely of the F = 7 — 6 hyperfine component, the B field estimates would be multiplied by a factor 
of 1.5. 

e S curve fits performed. 

Cross correlation method performed. Also, the line parameters (line center velocity, flux density) listed are 
for the strongest systemic feature in the spectrum. 

s Radial distance of maser feature from the central black hole. Positions derived from the line-of-sight velocities 
of the high- velocity features, assuming a Keplerian rotation curve and the standard values for JWbhi distance, 
systemic velocity of the galaxy and that the features lie on the midline. 

h Radial distance of maser feature from the central black hole. Position taken from Bragg ct al. (2000), which 
is based on the observed acceleration. 



the well-isolated features yield the same B field estimates 
and uncertainties. The systemic maser lines are stronger 
and more numerous, thereby improving the sensitivity to 
the B field, but have broader line widths than the high- 
velocity features and are heavily overlapped, decreasing 
the sensitivity. For the VLA data, we performed the 
cross correlation technique on those systemic maser lines 
that occupy the inner 90% of the bandwidth, in order 
to avoid passband problems at the edges of the spec- 
trum. We derived a value of B\\ = 29 ± 16 mG (Fig. EJ). 
Though this result is formally a 2 a detection, we do not 
count it as a significant detection (defined as a detec- 
tion with a significance level of 3 a), but rather report 
an upper limit of B|| < 50 mG (see Table |21 for com- 
putation). There are deviations near 493 km s _1 and 
at 499 km s _1 in the Stokes V spectrum, which are re- 
sponsible for the 2cr result. There is no clear detection 
of fractional linear polarization at a level of 0.3% (1 a) 
(Fig. El bottom panel). However, the conditions in the 
masing medium could give rise to internal Faraday ro- 
tation along the maser amplification path, leading to 
depolarizati on of a strong linearly polarized signa l (see, 
e.g.. 513.6 of lThompson. Moran. k Swenson 11200 ID . The 
requisite electron density for complete depolarization by 
this mechanism (i.e. differential Faraday rotation of n) 

is n e ss 7r (8.1 x 10 5 A 2 £?||-L) , where A is the observing 
frequency in meters, B\\ the B field along the line-of-sight 
in Gauss, L the pathlength over which depolarization is 
taking place in pc, and n e is the elect ron density in cm" 3 . 
For A = 1.35 cm, L = 0.003 pc ijMoran et alJ Il995j) . 
B\\ < 0.03 G, one obtains n e ~200 cm -3 . For an assumed 
average molecular density of 10 9 cm~ 3 , this corresponds 



to an ionization fraction of x e ~ 10 7 . These low val- 
ues are well within the range of the expected ion i zation 
frac tion. Models bvlNeufeld. Malonev. k Conger! <|1994D 
and iNeufeld k Malonevl 1)1995) , in which oblique X-ray 
illumination is heating the gas in the molecular accretion 
disk of NGC 4258 and other megamaser galaxies, give a 
theoretical upper limit of x e ~ 10~ 5 . Note that there are 
observations of strong linear polarization (up to 65%) i n 
galactic water masers ijGarav. Moran. fcH aschick 1989), 
which arise in star formation regions. 

Cross correlating the GBT spectra of the systemic fea- 
tures in LCP and RCP leads to a \-a upper limit of 
B|| <30 mG (Fig. Hand Table EJ. We chose the velocity 
range of the systemic maser lines such as to maximize 
the number of included features, while minimizing the 
extent of the baseline modulation seen in the renormal- 
ized Stokes V spectrum (bottom panel of Fig.[7J). Note 
that the different appearance of the systemic features in 
the two datasets are due to both time va riability and 
centr ipetal acceleration of ~10 km s^yr" 1 (Brag get alJ 
2000). 

In the case of the systemic features, Bu refers to the 
ordered radial component of the B field, assuming the 
features are threaded by a B field constant in magnitude 
and direction. They are located at a distanc e of R ~ 0.14 
pc (3 7, 000 Rs) from the central black hole ijBragg et alJ 
120001 as derived from the observed acceleration of the 
systemic features). The features are distributed over an 
azimuthal ang le in the disk of about 8°, or about 0.02 
pc. Recently. IHerrnstein et aD l|2005|) estimate that the 
systemic features are located in a bowl-shaped depression 
of the warped disk and occupy a region of length 0.05 pc 
along the LOS. 
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spectrum [Stokes V = (RCP-LCP)/2], 



Table [3 lists the details and summary of the individual 
Zeeman fits and their interpretation. 

4.3. Saturation Effects 

Some theoretical investigations have shown that the 
degree of observed circular water maser polarization de- 
pends also on the saturation state of the masers. Masers 
begin to saturate when R/T > 1, and reach full satura- 
tion when R/T > 100, where R is the rate of stimulated 
emission and T the decay rate, with r = 1 s^ 1 from the 
lifetime of involved IR transitions. 

The stimulated emission rate o f a maser is given b y 
R = A Su c 2 /(8 7T hv 3 ) x (0 6 /O s ) (IR.eid fc MoranllTflS^ . 
where A is the Einstein coefficient for spontaneous emis- 
sion (2x 10 -9 s _1 for the water transition), S v the specific 
flux density of the maser, and Of, and Q s are the beaming 
and subtended solid angle of the maser. For a cylindrical 
maser model, fib = (d/L) 2 , and Q s — (d/D) 2 , where d is 
the cross section diameter of the maser, L is the maser 
length, D is the distance to the maser source. Hence, we 
can write for D — 7.2 Mpc, 

S v \ ( L \ _j 



R = 50 



1 Jy 



10 



16 



cm 



(4) 



If the coherent velocity interval is 1 km s _1 , then the 
maximum coherent amplification length due to the line- 
of-sight gradient for a high-velocity maser feature at a 
radius of 0.2 pc is 0.2 pc x ^/2 km s _1 /u ro t or 3 x 10 16 
cm. Hence, i? m ax is 5 s _1 for the high-velocity masers 
with typical flux densities S v — 1 Jy. For the sys- 
temic masers, the beaming angle is known because of 



the clear velocity range of the systemic masers. It is 
8°, so that fib = 0.015 sr. In this case, we can write 



R 



0.8 S v (d/10 



16 



cm 



The maser spectrum is 



very crowded with the spacing between features of about 
1 km s _1 or larger. This means that d > 6 x 10 14 cm, and 
R is less than 500 for the systemic features with typical 
flux densities S v = 2 Jy. 

This means that the high-velocity masers are prob- 
ably not saturated. Furthermore, the linewidth of 
1.01 km s- 1 (FWHM) of the high-velocity feature at 
1269 km s _1 is inconsistent with a rebroadened, highly 
saturated maser that is a blend of the 3 strongest hyper- 
fine components, for the spectral line breadth at high sat- 
uration is expected to be FWHM = 1.5 km s~ 1 (N W92, 
lEmmering fc Watsonl Il99l iWatson & Wvldl I2003]) . In 
addition. iWatson fc Wvldl (|200C() favor R/T < 1 from 
numerical simulations of masers in a Keplerian disk with 
velocity gradients. The systemic features are probably 
slightly saturated, and might reach full saturation. 

Even if the masers are fully saturated, the upper limit 
on the m agnet ic field w ould be lower by a factor of max- 
imally 4 ijVOl Fig. 7; INW921 Fig. 7) - consequently, 
at any event, our reported magnetic field upper limits 
are truly upper limits. Since saturation increases the ob- 
served circular polarization, even in the event of strong 
saturation, the true magnetic field is smaller than the 
quoted upper limit. 



5. DISCUSSION 
5.1. The Total Magnetic Field 
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We obtained upper limits on the toroidal and on the ra- 
dial component of the magnetic field in the accretion disk 
of NGC 4258 by searching for Zeeman-induced circular 
polarization. Based on these results we did simulations 
to estimate the limit on the total magnetic field. Since 
the toroidal and the radial component, B tm and -B ra d, 
were probed at only slightly different radii, but different 
locations in the disk, we regard them as representative 
values for the disk and thus assume they can be used to 
characterize the magnetic field at a specific point, i.e., 
radius. 

We conducted a simulation by assigning random ori- 
entations to a range of total B field strengths, B to t = 
0— 1000 mG. For each value of Btot, 5000 random orien- 
tations were drawn, with uniform distributions in B tor , 
-B ra d, and B po \, the poloidal component. 

For each orientation, the value of the two components 
of interest were compared to our observational upper lim- 
its. Figure |H1 presents the results of the simulations: it 
shows the percentage of trials in which the toroidal and 
radial components of the intrinsic total B field could not 
be detected with our observational accuracy as a func- 
tion of the intrinsic total B field. We find that the simu- 
lated orientations for the B field would exceed measure- 
ments 50% of the time for a B to t value of 130 mG. Thus, 
we adopt the estimate of 130 mG as the upper limit on 
the total B field at a distance of 0.14 pc from the cen- 
tral SMBH. The Alfven speed, v A = B/y/Wp, would be 
2—20 km s _1 for the the allowable range of densities of 
molecular hydrogen, 10 s — 10 10 cm~ 3 , and for _Btot < 130 
mG. 

From simple hydrostatic equilibrium arguments where 
the magnetic pressure confines the masing gas, one could 
expect the following values for the B field: From B = 
y/8imn 2 kTk , the minimum value would be -B m ; n = 10 
mG (for the minimum temperature of Tk — 250 K 
and minimum molecular number density of nn 2 = 10 8 
cm~ 3 ) and the maximum value B max = 190 mG (for 
T k = 1000 K and n Ha = 10 10 cnr 3 ). LOS magnetic field 
measurements in starforming regions via water mascr 
polarization observation s lie within this range (FG89, 
ISarma et al.1 1200 ll I2002(K and polarization observations 
in envelopes of late-type starts indi cate super -thermal 
magnetic fields of up to - 0.1 - 3 G ifvo lfvol . 

An alternate way of estimating the total magnetic 
field is by considering predictions from simulations 
of magnetic fields in accretion disk. A series of 
magnetohydrodynamic (MHD) simulations show that 
the combined presence of magnetic field and out- 
wardly decreasing differential rotation can rapidly cre- 
ate MHD turbul ence via the magneto-r o tational in- 
stabi lity (MRI) llBalbus fc Hawleyl 119911: iStone et all 
IT99fl see also IBa.lhus Hawlevl 119981 IBalhusl 12003). 
This results in a greatly enhanced effective viscosity, 
which readily generates outward flow of a ngular momen- 
tum. Ha wlev. Gammie. fc Balbua l|199fifl performed lo- 
cal three-dimensional simulations of magnetized shearing 
box models, with an weak (p magn ctic < Pthcrmai) and 
zero-mean magnetic field as part of the initial condi- 
tions. Using a variety of initial field strengths and ge- 
ometries, they found that the large scale magnetic field 
is dominated by the toroidal component: B^ OI : B^ ad : 
-Bp ol « 10 : 5 : 1. If MHD turbulence is dominant, 



on e might expect a mag netic field configuration as found 
by IHawlev et al.1 l|199fifl . Hence, our measured upper 
limit on the radial component of the magnetic field then 
dictates an upper limit on the toroidal component of 
< V2 30 mG = 43 mG (GBT data), which is consistent 
with, and tighter than, the observationally determined 
upper limit of B tor < 90 mG at a slightly larger dis- 
tance (Table0). The MHD predicted total magnetic field 
would thus be -B to t = yl6/5 £? ra d S 50 mG at a radius 
of 0.14 pc, where we used the limit from the GBT data 
on the radial component, and £? t ot = \/16/10 -B t or HO 
mG at a radius of 0.2 pc, where we used the limit on the 
toroidal component. Since the radial limit is estimated 
closer to the black hole, we adopt Btot ^ 50 mG as the 
MHD predicted total magnetic field. 

There might be concerns that the magnetic field is not 
sufficiently coupled to the gas, as the ionization fraction 
prevalent in the molecular masing gas is probab l y very 
low. To probe this question, iMenou fc Quataertl 1)200 li ) 
conducted studies on ionization, MHD and gravitational 
instabilities in AGN accretion disks and found that the 
conditions in NGC 4258 's accretion disk still satisfy the 
requirements for large Reynolds numbers; hence, the neu- 
tral particles are expected to be well-coupled to the mag- 
netic field by collisions with the ions, and the MHD 
turbulence due to MRI in NGC 4258's accretion disk 
is expected to be similar to that in an ionized plasma. 
They also found that the thermal ionization instability 
is unlikely to be a source of efficient angular momentum 
transport and, thus, does not induce large-amplitude lu- 
minosity fluctuations in the disk. 

5.2. Effects of the Magnetic Field 

iReesI (119871): iBottqrff et all <|1997|) . and more recently, 
Bottorff & Fcrland (2000), have argued that magnetic 
fields of strength ~ 5 — 10 G confine the clouds that popu- 
late the broad-line regions (BLRs) in A GNs at di s tances 
R ~ 10 17 - 10 18 cm from the SMBH. IHo et all lfl997h 
report the detection of a broad Ha component in high- 
resolution, high S/N total- light spectra of NGC 4258 
(after decomposing the emission lines), with a FWHM 
of 1700 km s -1 . We note that the detected emission 
lines in polarized light (i.e., light scattered off material 
surrounding the central source) prob ably emanate from 
the narrow-line re gion, not the BLR (W ilkes et alJll995t 
IBarth et al.lll999T) . From the virial estimate, -Rblr, the 
position of the BLR, is approximately Mbh_G/v 2 , with 
v — /"fwhm, where the factor / depends on the details 
of the geometry, kinematics, and orientation of the BLR 
and is y / 3/2 for an isotr opic distribution with random 
inclinations l|Netzerlll99(t . Thus, i? B LR ~ 2.4 x 10 17 cm 
or 0.08 pc. This is smaller by a factor of 2 — 3 than the 
distance regime the magnetic field upper limits reported 
here are probing. In order to confine the BLRs at that 
distance, the B field strength would have to scale as r -4 
to reach the required high field strengths — since this 
scaling is unphysical, we conclude that magnetic pres- 
sure is unlikely to be the dominant force to con fine a 
BLR of radius 0.08 pc in NGC 4258 (see, however. iLaorl 
120031 for arguments on why NGC 4258 might not har- 
bor BLRs). However, since this system is an obscured 
AGN, the width of the broad Ha component could be 
larger, meaning the BLRs could be located closer to the 
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Fig. 8. — Result of simulations of the value of the total magnetic field (Btot) given our measured upper limits on the two components 
of the B field. The percentage of simulations where the total B field satisfies the upper limits on the toroidal and radial component and 
is therefore undetected with a sensitivity comparable to our measurements is plotted as a function of -Btot ■ The dashed line indicates the 
50% limit, i.e., in 50% of the trial runs the observational limits were met (see text for details). 



SMBH and thus high magnetic fields could be reached at 
those distances. We note that our upper limits reflect B 
fields in very dense masing gas and that the postulated B 
fields responsible for confining the BLR clouds would be 
presumably lower as they pervade a more rarefied inter- 
cloud medium (b e it relativis tic or accretion-driven wind 
or accretion flow; lReesHl987j) . 

Large scale poloidal fields may play an im- 
portant role in coUimating AGN jets (see, e.g. 
iBegelman. Blandford. fc Reesl Il98l . VLBI observa- 
tions of M87, which contains one of nearest large- 
scale jets, show that the jet is well-collimated at 
30 — 100 Schwarzschild radii, and are interpreted as 
suggesting an accretion disk threaded by a strong 
poloidal magnetic field as the launching mechanism 
i|Junor. Biretta. fc Livid I1999J) . The poloidal B field 
component in NGC 4258's accretion disk, which is 
observationally inaccessible, can be estimated to be 
-Bpoi ~ B lad /Vb < 14 mG (GBT data) at a radius 
0.14 pc (or 37,000 R s ) from MHD model constraints. 
Since for a steady-state thin accretion disk, which is in 
thermal equilibrium, the su rface density goes as r~ 3 / 4 
ijFrank. King. &: Raindl200^) . the upper limit on _B pol at 
r = 0.14 pc translates to B vo \ < 7 G at 10 Rs, for 
an a-thin disk that extends inwards to ~10 Rs- This is 
much less than the magnetic fields o f strength 10 2 — 10 4 
G req uired for coUimating radio jets ( Bland ford &: Pavnel 
1982). 



5.3. Mass Accretion Rate 



Relating the thermal pressure at the midplane of the 
accretion disk to the magnetic pressure, the mass ac- 
cretion rate can be estimated from the magnetic field 
stren gth. The thermal pre ssure on the midline is given b y 
(e.g-. lNeufeld fc Malonelll995llHe"rrnstein et alTll9 98b): 



1 



th 



pc s 



GMM 



37TV27T 



(5) 



where a is the Shakura-Sunyaev viscosity parameter, and 
M the mass of the black hole. 

Setting equation JSJ equal to f3B 2 /8ir, where the pa- 
rameter (3 is the ratio between the gas and the magnetic 
pressure f3 — Sirpc^/B 2 and is 1 for equipartition (note 
that (3 is sometimes defined as (3 — Ptot/Pm&g), we obtain 



M = 



3V2tt aPB^Csr* 



GM 



(G) 



If magnetic fields provide the kinematic viscosity, the 
parameters a and (3 are likely to be related according to 
the MHD-simulations prescrip tion a = C/(l + (3), w here 
C is usually set to 0.5 — 0.6 l)Naravan et al]ll99ll note 
their slightly different definitions of (3 and of the magnetic 
pressure, which are accounted for in the above equation). 
In that context, an equipartition magnetic field corre- 
sponds toa~ 0.3. Local three-dimensional MHD calcu- 
lations sugg est sub-thermal mag netic fields ((3 « 30—100, 
<(3> «60, lHawlev et al.l ll99fi). Recent global three- 
dimensional MHD simulations indica te a large mag netic 
field contribution, with (3 sa 1 — 5 (Hawlcy 2000J), be- 
cause the length-scales, where MRI can be effective, are 
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larger in the global simulations than in local ones. How- 
ever, they are conducted at radii very close to the SMBH 
(~ 1 — 10 Rs), owing to computational constraints, deal 
with torii at those distances (H/R ~ 1), and thus do not 
apply to our situation of a thin disk. 

Adopting a temperature of 400 K in the disk mid- 
plane, as suggested by analyses of observed water maser 
line profiles, equipartition between magnetic and ther- 
mal pressure, and a B field upper limit of 130 mG from 
the polarization analysis of the various maser lines and 
statistical considerations, we derive at the radius of r = 
0.14 pc a mass accretion rate of 



M< in- 3 - 7 °- 6 W 1 ) ( B tat > 
m ~ 1U l+(/3/l) V130 mG; 



,1/2 



(o.!4 pc) ( 



M BH 
3.9X10 7 M Q 



,400 KJ 



M yr- 



(7) 



If we a dopt the magnetic fi eld arrangement as pre- 
dicted by lHawlev et al.l l)1996|) with subthermal fields, 
such that B t 2 ot = 16B r 2 ad /5, and (5 w 60 (and thus 
a = 0.01), we obtain M < 10~ 4 2 M yr _:L at a radius of 
0.14 pc. We stress that this analysis is based on an upper 
limit on the average magnetic field acro ss the masing re- 
gion of ~ 7x 10 cm l|Moran et all! 995): if the magnetic 
fields are significantly tangled, then t he B field limit is 
unde restimated by a factor of a few ijWatson fe Wiebd 
2001, see also §OJ. 

The estimate in equation JJJ reflects an improvement 
by one order of magnitude from the previous estimate of 
the mass accretion rate in NGC 4258's disk (scaled to the 
presently estimated and measu red values of Mrh, Tk an d 
D) using the same methods (Herrnste in et al 



H, J-k auu 
L|[l998b). 



Our result is consistent with the mass accretion rate of 
10~ 4 a MQyr^ 1 currently estimated by J.R. Herrnstein et 
al. (2004, in preparation), which is based on considera- 
tions of the transition between mole cular and atomic gas 
using the X-ray irradiation model of[Ncufcld & Malonev 
( 1995) and incorporates new VLBI observations of the 
masers. 

The disk mass can be estimated by considering the de- 
viations of the maser velocities from pure Keplerian mo- 
tion (3 km s _1 ), which indicates that the mass pertaine d 
in the disk is less than about 1O 6 M ijMoran et all l9991. 
The associated average time-scale for accreting the whole 
disk is t = M^isk/M ~ 10 10 years, while the viscous time 
scale (i.e., the time scale for matter to diffuse through the 
disk from the masing region under the effects of viscous 
torques) is t v [ sc ~ R/vr with vr < 0.006a km s" 1 from 
the o bservations of the systemic features ijMoran et al.l 
1999, since Vr — av^,(H/R)), such that i v isc ~ 4 x 10 7 
years (for a = 0.3) and t v - lsc ~ 1 x 10 9 years (for 
a = 0.01). For comparison, the Eddington accretion 
rate for a central black hole mass of Mbh = 3.9 x 10 7 M Q 
is M E dd = 0.9 Moyr" 1 (e.g. IPetersoniri997l assuming 
a Shakura-Sunyaev disk), where the radiative efficiency 
parameter has the fiducial value of 0.1. Thus, our de- 
rived mass accretion rate is m = M/Mem ^ 10~ 37 for 
P = 1 and m < 10" 41 for f3 = 60. 

The above discussion relies on the assumption that 
the masers occupy the midplane of the accretion disk, 
for equilibrium considerations to be valid. There is 
no direct observational evidence for this claim, and 
in some models the masers inhabit a surface layer 



above a much colder and more massive accretion disk 
Ipesch. Wallin. fe Watsonl I1998|) or are uplifted by a 
disk-driven hydromagnetic wind ijKartie et alJ Il999j) 
(§ 15.41 gives a discussion of these alternate models). 
Nevertheless, the simple model of a warped Keplerian 
model and its cha racteristics fit remarkably well the ob- 



servational data (IHerrnstein, Greenhill. fc Morar] 



|Hjmmstgm|199^ 



1996, 



1999 



Humphre ys et al.l l2004). Any potential surface layer has 
to be very thin, H < 15 /iarcsec or 0.0005 pc (from 
the sp atial distribution of the masers in the syst ematic 
group, iMoran et all 1199.1 IHerrnstein et ai1l2005j) . For 
the systemic features, there is no direct way to mea- 
sure any potential displacement of the masers from the 
midline. We can attempt to provide a rough estimate 
by considering the radio continuum emission at 22 GHz 
from northern and southern lobes of the nuclear jet. If 
w e take the midpoint conn ecting the two lobes as shown 
in IHerrnstein et alJ l|1997|) as the location of the SMBH 
and compare it to the position of the dynamical center 
of mass as predicted by the maser kinematics, we esti- 
mate a difference of less than about 150 /iarcsec. If we 
take this limit as the offset between the maser layer and 
the midline (150 /iarcsec), i.e., 10 scale heights, we con- 
clude that the mass accretion rate upper limit (eq.[7]) is 
larger by 15 orders of magnitude, for a vertical pressure 
profile of P{z) oc (1 - Erf(z/\/2H)) 1 ^ 4 exp(-z 2 /2H 2 ), 
where Er f(x) i s the standard error function (see, e.g. 
IGoodmarJl200l . However, this estimate is most likely 
flawed, as the emission of the southern jet is probably 
attenuated by a layer of ionized gas in the disk via free- 
free or synchrotron absorption (He rrnstein et alJl!997j) 
and thus, the true centroid of the southern lobe may 
not be revealed in the 22 GHz maps. High-resolution 
submillimcter maps in the future era of ALMA could al- 
leviate the opacity problem and might answer the ques- 
tion fully. If the masers indeed arise on a surface layer 
and are not located on the midline, one would expect 
a discontinuity at the position of the black hole for an 
anti-symmetrically warped disk. The emission sites of 
the redshifted high-velocity features would be expected 
to lie on the top on one side, and the maser-skin of the 
blueshifted high-velocity features on the bottom on the 
other side. 

5.4. Comparisons to Maser Model Accretion Rates 

More than 50 AGNs with megamaser emission arc 
currently kno wn (L. Greenhill, p rivate communication, 
2004, see also lBraatz eTaTIISooa: iGreenhill et ai"l f2003'l. 
and several show clear signs of disk structure. NGC 4258 
is one of the best-studied cases and provides many ob- 
servational details, and numerous mechanisms have been 
proposed to explain the production of maser emission in 
its circumnuclear envir onment. 

Congerl l) 19941) and specifically 



Neufcld. Malonev. . 

INeufeld fc Malonevl 111995) propose irradiation of the 
warped circumnuclear accretion disk by X-rays from 
the cen tral AGN as the mechanis m to produce water 
masers. INeufeld fc Malonevl l|1995fl show that the outer 
edge of the maser distribution at r ~ 0.28 pc can be 
explained in terms of a phase transition from molecular 
to purely atomic gas, provided that the accretion rate is 



M w 10 



-4.1 



M yr 1 at a radius of 0.28 pc (see also 
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lMalonevl2002]) . Their mass accretion rate is in very good 
agreement with our upper limit. 

The inner radial cut-off of the maser annulus is postu- 
lated to be due the absence of the disk warp at the smaller 
radii and consequently, the absence of X-ray heating. 
Long-term monitoring observations of NGC 1068, an 
archetypical Seyfert II megamaser galaxy, show that the 
flux density of the nuclear masers are fl uctuating coher- 
ently on timescales of months to years ijGallimore et al.l 
120011) . These variations occur faster than the dynamical 
timescale would allow and can be invoked as response in 
reverberation to the central engine as the pump source, 
thus fitting the scenario of central X-ray heating (how- 
ever, it appears that the X-rays experience heavy ab- 
sorption and ca n only be detected in scattered light; 
IPier et all 11994^) . Similar monitoring campaigns of si- 
multaneous X-ray and maser emission observations are 
underway for NGC 4258 (L. Greenhill 2004, private com- 
munication) . 

As an alte rnative to the warped disk interp retation of 
megamasers, iKartie. Konigl. fc Elitzurl <[l999) propose a 
scenario where the masing gas is confined and uplifted 
from a flat accretion disk by a powerful, dusty, centrifu- 
gally driven MHD wind, which exposes the gas to X-ray 
irradiation and excites maser emission. This model pro- 
poses that the apparent warp seen in the maser disk of 
NGC 4258 is due to a non-coplanar distribution of mas- 
ing clouds and stems from their argument that the maser 
regions are clumpy, after considering radiative and kine- 
matic aspects and the time variability of the observed 
maser flux. The nearly perfect Keplerian shape of the ro- 
tation curve as traced by the maser features is explained 
qualitatively in terms of observational selection effects, 
but without fitting the wind-model to the observed spa- 
tial distribution of maser spots and their line-of-sight ve- 
locities. It is hard for this model to explain the anti- 
symmetric shape of the warp in the almost perfect Ke- 
plerian rotation curve of NGC 4258. Also, it is unclear 
how other observational trends can be explained with 
this model (e.g., the relative strengths of the two groups 
of high- velocity features, the apparent small height of 
the disk as traced by the masers). Their estimated line- 
of-sight B field of ~20 mG in the homogeneous outflow 
component at r ~0.2 pc is consistent with our upper 
limit of 90 mG at a similar position; furthermore, the 
B field in the masing cloud is expected to be smaller 
than the estimate in the outflow-wind, since the clouds 
are assumed to be diamagnetic (i.e., not permeated by 
an ambient magnetic field) in their model. Their derived 
mass accretion rate at r = 0.2 pc is M w 10 -11 M yr _1 
(their eq. [29]), assuming that 10% of the mass is car- 
ried away in the wind, and inconsistent with our value 
of M < 10~ 37 . Of course, a detailed comparison is not 
warranted here, since our accretion rate computation as- 
sumes that the masers lie in the accretion disk, the v ery 
opposite picture of the model bv IKartie et alJ l)1999t) . 

As an altern a tive t o the X-ray heating scenario, 
iMaoz fe McKed l|1998D envision spiral density waves 
sweeping through a self-gravitating accretion disk and 
triggering shock waves that heat the gas. This model 
is motivated in part to explain the consistently larger 
flux densities of the red-shifted over the blue-shifted 
high-velocity lines and the periodic spacing of the high- 



velocity features in radius and velocity. They predict a 
mass accretion rate of M « 10 -22 M Q yr _1 at the lo- 
cation of the masers (^0.2 pc) from energy dissipation 
arguments in the spiral density shocks. This is larger by 
almost two orders of magnitude than our present upper 
limit for the accretion rate at the same radius and thus, 
is not consistent with our results. Here, a direct compar- 
ison is warranted as the same disk geometry is assumed 
and prediction about the mass accretion ra t e con cerns 
the same location in the disk. iBragg et al.l l)2000[) find 
no observational evidence for the shock-model-predicted 
pattern in t he accelerations of the hi gh- velocity features. 
In addition, iHerrnstein et alJ l)1996t) show that the per- 
sistent asymmetry in the spectrum can be well- explained 
by a warped and X-ray irradiated disk. 

5.5. Comparisons to Accretion Models 

Provided accretion has been constant over the time 
span needed for material to reach the galactic nucleus 
from the maser disk, the estimate reported in this work 
can be used to probe the nature of the accretion respon- 
sible for the energy output of the central AGN. The vis- 
cous timescale is t v i sc ~ 4 x 10 7 years (for a — 0.3) 
or t v ; sc ~ 10 9 years (for a = 0.01; see § 15. 3|) and it is 
not certain whether the disk has been in steady state 
over the whole time range. For a jet whose evolution is 
closely linked to the evolution of the accretion disk, the 
dynamical life time of the jet and its relatively undis- 
turbed morpholog y gives an estimate of tj ~ 10 6 years 
(ICecil et al J 11995) as a time frame of steady accretion. 
The bolometric luminosity can be estimated from the X- 
ray luminosity, since Lx ~ O.lLboi- A number of X-ray 
observations have been conducted over the past several 
years and yield values for the absorbing column density 
of N H ~ (6 - 15) x 10 22 cm" 2 and for the 2-10 keV 
absorption corrected luminosity o f L x = (4 — 13) x 10 40 
erg s- 1 (for D = 7.2 Mpc) fsee iFruscione et al.ll200l 
Table 2 for a compilation and references therein). Thus, 
iboi ~ 10 42±1 erg s- 1 = 10- 3 ' 6±1 L E dd, where L Edd = 
1.3 x 10 38 (M B h/-/W©) erg s" 1 is the Eddington luminos- 
ity for this system (but see lWoo fc Urrvll2002t lLaorll2003l 
and for alternate estimates of the bolometric luminosity) . 

The highly sub-Eddington nature of the central source 
in NGC 4258 could be explained by an a disk that is ge- 
ometrically thin, optically thick and reaches to the cen- 
tral black hole and whose low mass accretion rate leads 
to the sub-Eddington luminosity of the central engine. If 
the observed nuclear jet as seen in the radio is powered 
by the accretion flow (and not by extraction of the spin 
energy of the SMBH), the minimum mass accretion rate 
has to be M ~ 1.2 x 10~ 3 eZ] &°/ 4 ( r/10) 3 4 Moyr" 1 
ijGammie. Naravan. fc Blandfordlll999j) . where j3 e is the 
ratio between the electron and the magnetic energy den- 
sity (assumed to be > 1 for jets) and T the Lorentz fac- 
tor of the relativistic electrons, which is > 5 (P. Kon- 
dratko, 2004, private communication, IHerrnstein et alJ 
1997). Thus, the accretion rate estimated in this paper 
would be sufficient to power the nuclear yet. We note 
that the standard thin disk model alone cannot gener- 
ate the large luminosity seen in X-ray portion of the 
SED, unless an X-ray emitti ng, hot, tenuous co r ona is 
present as, e.g. suggested bv lHaardt fc Maraschil l)1991l) 
for Seyfert galaxies. 
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Fig. 9. — Expected strength of the total B field as a function of mass accretion rate (eq. |7|) at a radius of 0.14 pc f or Tjy = 400 K, and 
for = 1 (solid curved line) and (3 = 60 (dashed curve d line). The accretion rates predicted by the jet-disk model of Falckc & Bicrmann 
( 1999) (JD AF), the thin dis k and the ADAF models of Gammi e et alJ 119991) are shown (dashed vertical lines), as well as the upper limit 
by Tflerrnstein et alj 11998b!) and that from the present experiment (horizontal lines). The vertically shaded region shows the constraint on 
the accretion rate assuming equipartition magnetic fields (i.e., Btot < 130 mG, /3 =1 and, thus a = 0.3) and the cross hatch shaded region 

depicts the constrai nt h aving incorporated the MHD predictions by Ha wlev et all ^996) (i.e., Btot = a/16/5 -B r ad J; 50 mG, /3 =60 and, 
thus a = 0.01, sec ij l5.1l for details). 



Alternatively, a number of authors have proposed 
that NGC 4258 is powered by an Advection- 
Dom ina ted Accretion Flow ( ADAF) llLasota et al.l 
Hflflfit IGammie et all Il999t ICharv et all 12000^ 
The ADAF model (e.g.. INaravan fc Yil 11991 see 
iNaravan. Mahadevan. fc QuataerTlT99ll ~and references 
therein for review) requires a geometrically thick, but 
optically thin disk, with a two-temperature plasma of 
hot ions (~ 10 11 K) and cooler electrons (~ 10 9 K). The 
high temperatures of the ions are maintained due to 
the inefficiency of collisional energy transfer (Coulomb 
interactions) at the expected high temperatures and low 
densities, i.e., accretion rates. Since the ADAF cannot 
cool radiatively, it "puffs" up (H/R ~ 1), and the ma- 
jority of the viscously dissipated energy gets transported 
by the ions into the black hole. The ADAF system emits 
at sub-Eddington luminosities due its low radiative effi- 
ciency of the energy-carrying ions and has been put forth 
in explaining low-luminosity AGNs (such as NGC 4258), 
which might occupy a significant fr acti on of the local 
galaxy population l)Ho et alJ E?95j). IGammie et al.1 
( 1999) fit an ADAF model to the spectral energy 
distribution (SED) of NGC 4258 wh ile incorporating 
the stringent 22 GHz upper limit by iHerrns tein et all 
(|1998aj) and the IR data bv ICharv fc Becklinl lj!997). 
which constrain the transition radius between ADAF 
and thin disk to lie at 5 - 50R S = (10 - 1 00) GM/c 2 . 
Their derived mass accretion rates for slightly different 
ADAF models cluster around M = 10~ 2 M Q yr~ 1 . This 



ADAF accretion rate is not consistent with our reported 
accretion rate. 

Figure [§] shows a comparison of the different M val- 
ues including our upper limit on the total B field 
from our measurements in conjunction with statisti- 
cal considerations and from MHD predictions. Detec- 
tion of a relativistically broadened Ka emission line 
at 6.4 keV would conclusively demonstrate the pres- 
ence of cold material close to the last stable or- 
bit in NGC 4258 and rule out the standard ADAF 
model. Detection claims have been made of a narrow 
iron component ijRevnolds. Nowa k. fc Malonevl 120001 ) 
that is apparently tim e-variable ijPietsch & Read l2002t 
iFruscione et alJ 120051) . but there has been no detection 
of a broad emission component. 

Nonetheless, since both ADAF and thin disk models 
predict a spectral behavior of F v oc v 1 ^ over the IR por- 
tion of the SED (from the emission of the thin disk), 
they are in conflict with t he most rece n t high -resolution 
mid-IR data obtained bv ICharv et al.1 l|2000D that indi- 
cate a non-inverted and a much steeper power-law (F v oc 
v~ 1A ). A jet-dominated-accretion-flow, which has been 
proposed as an a lternative to both the ADAF and the 
thin disk model l|Falcke fc BiermanrJ 119991 I Yuan et alJ 
2002), is able to accommodate the large IR flux and its 
spectral behavior. The compact jet that is evoked to give 
rise to the whole extent of the SED would be fueled by 
an outer ADAF and a further-out lying thin accretion 
disk. The mass accretion rate that is needed to power 
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the jet appears to be low: iFalcke fe Biermanr] l)1999fl 
esti mated an accretio n rate of the order 10 -4 M©yr _1 
and lYuan et alJ l)2002j) showed that jet mass loss rates of 
~ (0.7-2) x 10~ 4 M0yr _1 generate viable JDAF spectra 
in agreement with data over a broad range of energies. 
The exact relationship between jet l oss and disk accre- 
tion rate is uncertain, nevertheless, I Yuan et alJ l)2002f) 
show that low accretion rates (M ~ io~( 4 ~ 5 ) M©yr -1 ) 
are required in order to have negligible emission from the 
disk and the ADAF. These values are are consistent with 
our estimated accretion rate (Fig. |j|J). 

There exist a number of other accretion scenar- 
ios that belong to the larger class of radiatively in- 
efficient accretion flows (which include ADAFs and 
JDAFs) and that could explain the accretion pro- 
cesses operating in the central engines of low- 
luminosity AG Ns: Convection Dominated Accretion 
Flow s (CDAFs ;|Naravan, Igum enshchev. fc Abramowicd 
1200(11 iQuataert fc Gruzinovl 12000]) and Adiabatic In- 
flow Outflow Solutions fADIOS: iBlandford fc Begelmanl 
1999). As these models have not been fit to the spectrum 
of NGC 4258 in detail, they will not be discussed here 
at length. Qualitatively, the CDAF model is expected to 
be mass-starved, since the inward transport of angular 
momentum by convection chokes accretion. The ADIOS 
model stipulates a powerful wind, which carries away en- 
ergy and angular momentum from the accreted gas; thus 
it would require a large accretion rate at large distances 
from the SMBH to compensate for the mass loss occur- 
ring further in, leading to a low mass accretion near the 
event horizon. 

6. SUMMARY AND CONCLUSIONS 

In this paper, we have presented total power and polar- 
ization spectra of the nuclear maser in NGC 4258 taken 
with the VLA and the GBT. These spectra constitute 
the most sensitive data with the most complete simulta- 
neous velocity coverage of maser emission in this system. 
Based on the analysis of the polarimetric observations we 
conclude the following: 

1) There is no circular polarization indicative of 
Zeeman-induced splitting of the maser lines of water va- 
por, a non-paramagnetic molecule, at the level of 0.1 — 
0.3%. This leads to a 1 a upper limit on the toroidal 
component of the magnetic field of 90 mG (at 0.25 pc, at 
the position of the redshiftcd high- velocity features) and 
on the radial component of 30 mG (at 0.14 pc, at the po- 
sition of the systemic features) in the maser-delineated 
accretion disk of NGC 4258. No linear polarization at a 
level of ~0.1% is detected. If the masers are saturated, 
the B field upper limits may be lower than reported here. 

2) The magnetic field upper limits are less than the 
several-Gauss-fields required to magnetically confine the 
broad line regions clouds, whose emission is seen in high- 
resolution, high signal-to-noise optical data in NGC 4258 
l)Ho et al.ll997|) . provided that the clouds are close to the 
position of the maser features (as suggested by the virial 
estimate). 

3) Assuming equipartition between magnetic and ther- 
mal energy in the accretion disk (and thus a = 0.3) 
and that the masers are situated on the midline, we 
derive a mass accretion rate of < 10 _3 7 M©yr _1 for 
a total magnetic field of less than 130 mG. Alterna- 
tively, adopting the B field configuration as predicted 



by MHD simulations of lHawlev et alJ l)1996|) . we obtain 
M < 10 _4 2 M©yr _1 . Both upper limits on the accretion 
rate are in agreement with estimations from the X-ra y 
irradiation maser models bv lNeufeld fc Malonevl fl995), 
but do no t agree with the high M from the spiral density 
model by Maoz & McKee (1998) predicted at the same 
radius. 

4) Provided the accretion disk has been in steady state 
and the masers reside on the midline, the low accretion 
rate reported here favors mass-starved models as candi- 
dates for explaining the extremely low-luminosity nature 
of NGC 4258. These models include the thin-disk model 
and the jet-disk model. 

Detection of the magnetic field in the accretion disk of 
NGC 4258 could be achieved with more sensitive data. It 
would be especially important to obtain detections in the 
numerous high- velocity features in order to map the mag- 
netic field structure as a function of radius in an AGN 
accretion disk. MHD theories could be tested, as they 
predict random, not ordered, magnetic fields at different 
radii. Detection of the radial component of the mag- 
netic field via the systemic features might be hampered 
by saturation effects. 
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APPENDIX 

COMPARISON TO THERMAL ZEEMAN EFFECT TERMINOLOGY IN RADIO ASTRONOMY 

Another instance of fractional Zeeman splitting can be found in thermal OH or H I emission or absorption obser- 
vations. Even though the Zeeman splitting coefficients are ~ 10 3 times larger for those paramagnetic molecules than 
for water, the regions in which the emission originate are of low density, such that the B fields are in the range of /iG, 
hence producing small Zeeman splittings. Numerous thermal Zeeman effect radio experiments thus fit the following 
expression to the Stokes V spectrum (see ITroland fc Heileslll982t ISault et aLlll990l for good reviews) : 

V(v) = Cl ^P- + c 2 I(v). (Al) 

Here, C2 denotes the instrumental polarization term due to small calibration errors between LCP and RCP. The factor 
ci is a measure of the Zee man sp litting and corresponds to c\ = Ap-F'-Bcos^/v^- This notation is equivalent to the 
terminology introduced by FG89 (eq. [Q), as shown in the following paragraph. 

We assume that the Stokes / profile can be expressed as a simple Gaussian line: Stokes I — I e~ v l 2a (for simplicity, 
the center velocity is taken to be at km s" 1 ). For the small Zeeman splitting case , the Stokes V pro file is proportional 
to the Stokes I profile, or more specifically (;Goodmaril ll989t ICrutcher et al]ll993tlHerrnsteinlll997|) . 

^ 1 ~dv~ 

= I a ^veS, (A2) 

Thus, the S curve is an antisymmetric curve and has extrema at ±<r (not at ±a/y/2 as claimed bv lVlemmings et alJ 
120011 |2002|) with values of ±(I oCl /a) e" 1 / 2 . 



Hence, 



Kiax Kiax — 1/2 



^=EL Bcos0 ±21. (A3) 
V2 a 



Since the FWHM of a Gaussian profile is defined as FWHM = 2oV2 ln2 = Avl, we can substitute a and simplify: 



Vmax 
-^max 



Thus, we have arrived at equation gTTp . 

CROSS CORRELATION METHOD 

When the magnetic field is constant across the spectrum we can estimate the Zeeman splitting across an ensemble 
of lines by cross correlating the LCP and RCP spectra without forming the Stokes V spectrum. Here we show that 
the sensitivity of this method is equivalent to that of the S Curve method. 

Determining the net velocity shift between the LCP and the RCP spectra can be achieved by simply finding the 
maximum in the cross correlation function. Since the Zeeman splitting is very small for the expected sub-Gauss B 
fields, the maximum of the lag function will lie between lags of and +1 channels or lags of and —1 channels. 
Assuming a parabolic shape for the cross correlation function at maximum, we can write the following 3 equations: 

CC+ =a{A-q) 2 + b 
CC = a (q) 2 + b 

CC- = a (-A - q) 2 + b, (Bl) 

where A denotes the channel spacing (in km s _1 ) and q the velocity at which the cross correlation function peaks. 
Here, CC+ is the value of the cross correlation function at a lag channel of +1, i.e., 

N 

CC + = £>y i+1 , (B2) 

i=l 

with Xj denoting the value of the LCP spectrum in the \ th channel, y; denoting the value of the RCP spectrum in 
the i th channel and N the number of channels. Analogously, CCq corresponds to the value of the cross correlation 



^4f-F' n a 2V2 ln2 _ 1/2 

— =— BcosQ — e ' 

V2 Av L 

A F -F^COsfl 



18 



Modjaz et al. 



function at a lag channel of 0, and CC_ at a lag channel of —1. Solving the 3 equations with 3 unknowns, the net 
velocity shift q can be obtained with 

= A CC + - CC. 
q 2 CC- + CC+ - 2 CC ' K ' 

We can obtain the conversion coe fficient f rom velocity shift to magnetic field from the synthetic spectra (similarly 
to the 5*-curve method calibration bv lFG89f) . with 

B " = Trib' (B4) 

where Af-f' — 0.013 km s" 1 G -1 for the F = 7 — 6 transition. However, we adopt the value derived from radiative 
transfer calculation bv lNW92l of A F - F , = 0.020 km s" 1 G~\ appropriate for the merging of the 3 strongest hypcrfinc 
components for masers. 

The Uncertainty in Bn from the Cross Correlation Method 

In equation l|B2|) . each channel value in the spectrum is a measured variable. We derive the uncertainty in the 
estimate of q by using the standard error propagation equation: 

2 \p / dq , 2 2 , / d 1 n2 2 , 9 9q dq 



Since the LCP and RCP receivers are independent, the noise in the LCP and RCP spectra are not correlated, and 
e covariance 
Now, define 



the covariance term a xy can be set to zero 



f = cc+- cc- = Eti( x i yi+i - x i yi-i) ( B 6) 

g =CC- + CC+ -2CC = Eti(xi Yi-i + xi y i+ i - 2 Xi yi ). (B7) 



Thus, 



9f 3g r 

dq_ = A a^g ~ gxj 1 = A (y j+ i - yj_i) g - (yj_i + y j+ i - 2 yQ f 
<9xj 2 g 2 " 2 g 2 

A (xj_i - x j+ i) g - (xj_i + x j+ i - 2 Xj) f 
5 Yj "2 g 2 



(B8) 
(B9) 



Finally, 



j=i J 
A 2 w 

=tiE - ^j- 1 ) s - (yj- 1 + y^ 1 - 2 ^) f ] 2 ^ + ( Bn ) 

3 3=1 

[( Xj _ x - x j+1 ) g - ( Xj _! + x j+1 - 2 Xi ) f ] 2 ^. (B12) 

Note that a x and are the rms noise in the LCP and RCP spectra, respectively, and that they were derived using 
line-free channels. 
Thus, from equation (|B4|I . 

- VIST' <B13) 



Monte Carlo Simulations 

We performed Monte Carlo simulations to understand the proper implementation of both the S'-curve and the cross 
correlation method, and to ensure that they have the same sensitivity. During these simulations we generated synthetic 
Zeeman S spectra as produced by a strong B field. In order to simulate data, we injected repeatedly typical noise 
levels to the synthetic spectra and sampled them at typical channel spacings (here 0.164 km s" 1 , see §|3J|. We subjected 
the ensemble of spectra to the same analysis procedures as the real data set, measuring B||. The uncertainty for the 
individual B\\ field "measurement", as u , was derived from a formal minimum x 2 fit for the iS-curve method and from 
error propa gation considerations for the cross correlation method. 

As Figurc lBlOl shows. both methods can recover the input magnetic field strength and both have similar uncertainties. 
The uncertainty in the B\\ field estimation is expressed in two ways: a) the width of the histogram of the Gaussian 
distributed Bn field measurements (top panel, Fig. IB10|) displays the intrinsic uncertainty in performing Zeeman 
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measurements on a dataset with the characteristic noise and data qualities; b) the value for the center of the distribution 
of cr_B||, the -B|| field uncertainties (bottom panel) as reported by the two Zeeman methods. Since the values for a) 
and b) coincide for both methods, we are confident that i) the newly developed cross correlation technique is properly 
measuring the Zeeman splitting and ii) our estimation of the uncertainties for both methods are accurate. With these 
simulations, we are also able to get a handle on the uncertainty of the uncertainty in the B\\ field estimate, <Vb , 

indicated by the width of the histogram for the measured Bn field uncertainties (bottom panel). Again, the values for 
a„ B coincide for both methods, as expected. 




500 



400 



300 



200 



100 



Center of Distribution = 
Gaussian width 



43.51 mG 
= 6.18 mG 




20 30 40 50 60 70 

Value of (j g from cross correlation method [mG] 



500 



400 



300 



200 



100 



Center of Distribution = 
Gaussian width 




42.08 
■ 5.46 



20 30 40 50 60 70 

Value of <J B from S curve method [mG] 



Fig. BIO. — Results of 10,000 Monte Carlo simulations of the 5-curve and the cross correlation method. They were performed for 
synthetic data that consist of a Gaussian line with an amplitude of 1212 mjy and a line with of 0.79 km s -1 (FWHM) plus an rms noise 
of 2.9 mjy (or a resolution of 0.16 km s _1 ). The histograms show the distribution of "measured" magnetic fields B and their associated 
uncertainties cr^ . Superimposed are Gaussian fits to the histograms and their parameters arc displayed. Both methods arc able to recover 
the applied B field over the statistical range of simulations and both have similar uncer tainties. Note that the parameters of the simulated 
data shown here are those of the GBT high- velocity feature at 1303 km s _1 (see S I4.1I . 
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